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Measurement of triple gauge WW7 couplings at 
LEP2 using photonic events. 



A study of events with photons and missing energy has been performed with the 
data sample obtained with the ALEPH detector at centre-of-mass energies from 161 
to 184 GeV, corresponding to a total integrated luminosity of about 80 pb _1 . The 
measured distributions are in agreement with Standard Model predictions, leading 
to constraints on WW7 gauge coupling parameters Ak 7 and A 7 . The results from 
the fit to the cross sections and to the energy and angular distributions of the 
photons are: 
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1 Introduction 

In e + e~ collisions at LEP 2 energies, the trilinear WW7 and WWZ couplings can be 
probed with direct W-pair (e + e~ — > W + W~), single W (e + e~ — > Wei/) production or 
with photon production (e + e~ — > 1/1/7(7) ) [1,2], In the WW channel a minimal set of 
five independent parameters is necessary to describe the Z and 7 couplings to the W, 
assuming C and CP conservation. Usually a model-dependence is introduced to reduce 
this set to at most three parameters (e.g. the model with the parameters aw, ®-W4>, a B<t> 
||). Although the photonic channel is less sensitive to the couplings than the W pair 
and single W channels |4j], it can resolve sign ambiguities and is therefore complementary. 
Constraints on the WW7 vertex have also been obtained at the Tevatron [H within a 
slightly different theoretical framework. 

The purpose of this letter is to set constraints on the WW7 trilinear couplings with 
a study of photonic events, using data collected by ALEPH at centre-of-mass energies 
ranging from 161 to 184 GeV and corresponding to a total integrated luminosity of about 
80 pb" 1 . 

In the Standard Model, three processes contribute at tree level to the 1/F7 final state 
corresponding to the five diagrams shown in Figure 1. 
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Figure 1: Feynman diagrams for e + e — > 1/1/7. Only the process C is sensitive to the 
WW7 couplings. 

The WW7 vertex is present only in the last diagram, which contributes about 0.3% to 
the total Standard Model e + e~ — > vvy cross section, but which also leads to characteristic 
energy and angular distributions of the final state photons. A measurement of the total 
cross section, supplemented by a fit to these distributions, is therefore sensitive to the 
presence of an anomalous WW7 coupling. 

This vertex can be described by three C and P conserving parameters, gl, k 7 and 
A 7 , related to the following W boson properties: 



charge Q w = egl 



7 



magnetic dipole moment fx w = {gl + n y + A 7 ) 



In the Standard Model, these three parameters are equal to 1, 1 and 0, respectively, 
and their deviations from these values are parameterized as "anomalous couplings" Ag 7 , 
Ak 7 and A 7 . Here, the electric charge of the W boson is assumed to be equal to that of 
the electron, thus fixing — 1, while no further assumptions are made on k 7 and A 7 . 
The matrix element for the e + e~ — > 1/1/7(7) final state is a linear function of Ak 7 and 
A 7 . Its implementation in the KORALZ Monte Carlo program |J, including initial state 
radiation of additional photons, is used throughout this analysis. 

This letter is organized as follows. In Section 2, the aspects of the ALEPH detector 
relevant to this analysis are described. The event selection is presented in Section 3, and 
the fit of the data to the presence of anomalous WW7 couplings is discussed in Section 
4. Section 5 gives the fitted values and the resulting constraints on Ak 7 and A 7 ; the 
systematic uncertainties are discussed in Section 6. 

2 The ALEPH detector 

The ALEPH detector and its performance are described in detail in 0, U- Here only a 
brief description of the properties relevant to the present analysis is given. 

The central part is dedicated to the detection of charged particles. From the interaction 
point outwards, the trajectory of a charged particle is measured by a two-layer silicon 
strip vertex detector, a cylindrical drift chamber and a large time projection chamber 
(TPC). The three tracking detectors are immersed in a 1.5 T axial field provided by a 
superconducting solenoidal coil. 

Photons are identified in the electromagnetic calorimeter (ECAL), situated between 
the TPC and the coil. It is a lead-proportional-wire sampling calorimeter segmented in 
0.9° x 0.9° towers read out in three sections in depth. It has a total thickness of 22 radiation 
lengths and yields an energy resolution of 5E/E = O.lS/yE + 0.009 (E in GeV). Two 
independent readouts of the energy are implemented respectively on the cathode pads and 
on the anode wires of the ECAL. At low polar angles, the ECAL is supplemented by two 
calorimeters, LCAL and SiCAL, principally used to measure the integrated luminosity 
collected by the experiment, but used also here for vetoeing purposes. 

The iron return yoke is equipped with 23 layers of streamer tubes and forms the 
hadron calorimeter (HCAL), seven interaction lengths thick; it provides a relative energy 
resolution of charged and neutral hadrons of 0.85 /\/~E. Muons are identified using hits 
in the HCAL and the muon chambers; the latter are composed of two layers of streamer 
tubes outside the HCAL. 

The information from the tracking detectors and the calorimeters are combined in an 
energy flow algorithm ||. For each event, the algorithm provides a set of charged and 
neutral reconstructed particles, called energy flow objects, used in the analysis. 

3 Event samples and selection 

The data were collected with the ALEPH detector at LEP at several centre-of-mass 
energies between 161 and 172 GeV in 1996, and between 181 and 184 GeV in 1997. 
The corresponding integrated luminosities are given in Table [I]. 



Table 1: Data samples. 



Energy 
(GeV) 


Luminosity 
(pb- 1 ) " 


N events 
Data Expected 


161 


11.0 


32 31.8 


172 


10.7 


27 32.2 


183 


58.1 


148 145.8 


Total 


79.8 


207 209.8 



3.1 Selections and cuts 

Photon candidates are defined as described in ||. Only events with no reconstructed 
charged particle tracks and at least one photon with an energy E 1 > 0.1y/s are considered; 
the trigger efficiency for such events is almost 100%. At most one hit is accepted in the 
muon chambers, to eliminate beam-related and cosmic ray muons. The loss of signal 
events with noisy muon chambers was estimated from events triggered at random beam 
crossings to be 3%. The timing of the energy deposition in the ECAL is checked to be 
consistent with the beam crossing time. 

All events with at least 0.5 GeV detected below 14° from the beam axis are rejected, 
in order to remove radiative Bhabha events. The efficiency correction factor associated 
with this cut was estimated from events triggered at random beam crossings to be 3.5%. 

The consistency between the energy measured from the ECAL pads and from the 
ECAL wires is checked. In case of leakage out of the ECAL, a localized energy deposit in 
the HCAL, E had , associated to an ECAL cluster is added to E 7 , after correcting for the e/ ti 
ratio; only events with E^/E^ < 10% are kept. To reduce the remaining background, 
all but 2.5 GeV of the total energy is required to come from photon candidates. 

At least one photon candidate is required to fulfil the conditions 9 1 > 20° and 
Pr-y/Ebeam > 0.1. For multiplioton candidates, the additional photons are considered only 
if their energy exceeds 0.05^. The overall missing transverse momentum is required 
to be greater than 12 GeV/ c. The last cut removes the remaining Bhabha events with 
radiation at large angle. 

Table |I] shows the data samples used in this analysis. The numbers of selected 
events agree with the numbers expected from the SM cross sections determined with 
the KORALZ Monte Carlo. The cross section measured from the data at 183 GeV, with 
the present analysis and within the global kinematic cuts, is 3.45±0.30 pb, to be compared 
to the SM prediction of 3.40 ± 0.02 pb. 

3.2 Monte Carlo simulation with KORALZ 

The simulation uses a modified version of the KORALZ program, which includes the SM 
expectation (with electroweak corrections) as well as QED radiative corrections, and the 
contributions of anomalous coupling amplitudes with exact matrix element calculations 
||. The overall higher-order QED correction factor is around 1.4, but depends on the 
centre-of-mass energy. More than ten thousand simulated events are used for each energy. 
To obtain a description of the anomalous couplings in the simulation, each event 



is assigned a weight, which is a function of Ak 7 and A 7 . This method provides the 
smallest uncertainty, as the statistical error corresponds only to the differences between 
the distributions produced from the Standard Model and those from anomalous matrix 
elements. 

As the matrix element is linear in Ak 7 and A 7 the cross section and the differential 
distributions are bilinear forms of Ak 7 and A 7 . For each event it is thus sufficient to store 
weights for only six configurations in the (Ak 7 , A 7 ) plane, in order to compute any cross 
section or kinematic variable as a function of Ak 7 and A 7 . The leading order amplitudes 
including these anomalous couplings are folded with higher order QED effects, following 
the procedure discussed in JTTJ. 

The Standard Model predicts that the cross section for the radiative return to the 
Z resonance decreases when the centre-of-mass energy increases, while the opposite is 
true for the W exchange. In case of anomalous contributions, the sensitivity of the cross 
section increases almost quadratically above 161 GeV. 

The kinematic cuts have been chosen to optimize the sensitivity to the anomalous 
couplings Ak 7 and A 7 . Figure || gives the statistical sensitivity, defined as the anomalous 
contribution divided by the statistical error on the SM expectation, as a function of the 
scaled energy variable xe = E 1 jEy )eam . The minimum of the sensitivity occurs around the 
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Figure 2: Statistical sensitivity of this analysis as a function of the scaled energy xe at 
183 GeV centre-of-mass energy: a) for Ak 7 , b) for A 7 . The sensitivity is defined as the 
anomalous contribution divided by the statistical error on the SM expectation. The solid 
and dashed histograms correspond to parameter values of —10 and +10, respectively. The 
single photon radiative return to the Z corresponds to xe = 0.75. 
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Figure 3: Inclusive distribution of a) the scaled energy xe and b) the absolute value of 
the cosine of the polar angle of the photons, after all selections, for data and Monte Carlo 
at 183 GeV. 

position of the Z return peak (xe = 0.75). An important observation is that for Ak 7 > 
the differential cross section decreases almost linearly as a function of Ak 7 for events with 
xe < 0.75, and increases quadratically above. 



4 Likelihood fit 

In addition to the observed number of events, two kinematic variables of the photon are 
used in the fit: the photon polar angle 6 1 and the scaled energy xe- Figure |3| shows the 
distribution of xe and | cos# 7 | for data, compared to the Standard Model predictions for 
y/s = 183 GeV. 

Two kinematic regions have been chosen for the fit, excluding the region of the Z peak 
return where the sensitivity to the anomalous couplings is minimal. Only photons with 
| cos# 7 | < 0.9 are used for the fit to the shape of the distributions. 

Defining E'z = [s — m|)/2 v /i, the two kinematic regions are the following: 

• Region 1, low energy photons with E 7 < Ez — 3Fz 

The contribution from higher order radiative corrections is described by an almost 
constant term obtained from the Monte Carlo simulation. The scaled variable xg is 



found to be as discriminant as the angular variable in the fit. Both are used for the 
Ak 7 fit, whereas A 7 is determined only from the total cross section. The sensitivity 
to the A 7 parameter in this kinematic region is very low. 

• Region 2, high energy photons with E 1 > £!? + 0.5 GeV 

In this region, the higher order radiative corrections decrease the number of expected 
events by 30%. The scaled energy xe is more discriminant than the angular variable, 
both variables being used in the fit of Ak 7 and A 7 . It can be observed (Figure 0) 
that the sensitivity to A 7 with xe is similar to that of Ak 7 . 

Limits for anomalous coupling parameters have been derived from the generalized 
likelihood expression: 

logL = log 1 thJ + log ( th) (2)| + E logP[ 1} + E logPP, 

obs' obs ' 

where P- , are the probability density functions of observing event i with a given 
value of xe and 9 1 in region 1 and 2 respectively, and N^j and are the expected 
number of events in each region, including background. This likelihood formula contains 
two parts: the first one concerning the number of observed events, the second one being 
related to differential distributions for each kinematic region. The number of events used 
in the fit and those expected from the SM are given in Table 2. 

The acceptance convoluted with the experimental resolution leads to correction factors 
to the cross sections of 1.10 for the first kinematic region and 0.7 for the second; these 
correction factors are constant (within ±2%) in each region as Ak 7 or A 7 vary. 

The studies made with the KORALZ Monte Carlo show that the cross sections and 
distribution shapes vary differently in the two kinematic regions. For low energy photons 
the anomalous effects result from the interference term between the SM amplitude and the 
anomalous amplitude; the resulting variation is monotonic and linear for Ak 7 (A 7 ) > 
and Ak 7 (A 7 ) < and only one solution is expected for the A/t 7 and A 7 fit. For the 
high energy photons, the variations are quadratic (due to a quadratic contribution of the 
anomalous amplitude) and one or two solutions are expected; the case of one solution 
corresponds to A/t 7 = or A 7 = 0. This behaviour, important in the error determination, 
is discussed later when the error calibration procedure is presented. 



Table 2: Number of events (N Events) entering the fit in the two kinematic regions. The 
number of expected events is estimated from the KORALZ cross sections, corrected for 
acceptance. 



Kinematic 
region 


N Events, 
Cross section fit 
Data Expected 


N Events, 

(X E ,0) fit 

Data Expected 


Region 1 
Region 2 


93 101.0 
30 32.8 


60 67.4 
23 25.6 




Figure 4: Likelihood curves for the fit of Ak 7 at A 7 = for the contribution of the cross 
section term (solid curve) and the shape term in and 9 (dashed curve). 

5 Results 

The likelihood functions are calculated globally for the cross section and on an event- 
by-event basis for the energy and angular distributions. Figure |4] displays the variations 
of the log-likelihood (— AlogL) corresponding to the fit of Ak 7 at A 7 = 0, for the cross 
section and the distribution contributions. At present energies, the contributions of the 
cross section and of the shape variation terms are equally important for the fit of Ak 7 . 
The result for A 7 is dominated by the sensitivity to the shape in Region 2. 

Figure |5| shows the (—AlogL) functions for Ak 7 fitted at A 7 = 0, and for A 7 fitted at 
Ak 7 = when the two contributions are merged. The results are: 

Ak 7 = 0.05^J;j;Q(stat) assuming A 7 = 
A 7 = — 0.05^};45(stat) assuming Ak 7 = 

where the errors correspond to an increase of — logL by 0.5. The lower precision for A 7 is 
expected since the exchanged W's are at a rather low momentum scale and the A 7 term 
in the Lagrangian contains high powers of the W momentum. 

The 95% CL. limits derived from the one parameter fits are : 

—2.1 < Ak 7 < 2.2 assuming A 7 = 
—3.0 < A 7 < 3.1 assuming Ak 7 = 0. 

The validity of these 95% CL. limits have been checked using 100 Monte Carlo samples 
corresponding to the data luminosity, the analysis procedure described for the data being 
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Figure 5: Likelihood curves for the fit of A 7 at Ak 7 = (solid curve) and A/t 7 at A 7 = 
(dashed curve) for the sum of the cross section and distribution shape terms. 

applied to each Monte Carlo sample. This study indicates that these errors are consistent 
with the frequentist interpretation, within 10% of their values, and do not benefit from 
favourable statistical fluctuations. 

Figure |B| shows the 68% and 95% confidence level contours in the ( Ak 7 , A 7 ) plane from 
a two-parameter fit. Although the two parameters are not independent, the confidence 
level contours are symmetric. This comes from the fact that the results are very close 
to 0, so that only one minimum is found. If the results were far away from the SM 
prediction, there could be several local minima, around which the two parameters would 
be correlated. 

6 Systematic uncertainties 

The contributions to the systematic uncertainty on the determination of Ak 7 are 
summarized in Table [3|. The total systematic uncertainty is much smaller than the 
statistical one. 

• The acceptance corrections were tested with different cuts in xe and 9. This led to 
an uncertainty on the fit results as shown in Table 3. 

• The main contribution to the systematic error in the present study comes from the 
energy calibration of high energy photons, which has been checked to be 1% with a 
large sample of e + e~ — > 77 events. 
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Figure 6: 68% and 95% confidence level contours in the Ak 7 , A 7 plane. 



The possible contributions to the e + e~ — > 7(7) + X channel, other than X = uV, 
may come from radiative Bhabhas or e + e~ — > 77(7) events. All such events in the 
Monte Carlo sample are eliminated by the angular and energy cuts. 

The KORALZ simulation of higher order effects gives a correction of about +100% 
for the SM cross section in Region 1, and about —30% in Region 2. 

A theoretical estimate of the error on these correction factors is about 5%. However, 
only comparisons with complete calculations from the exact matrix elements (not 
present in KORALZ) for the two and three hard bremsstrahlung photons would 
allow a satisfactory estimation of this uncertainty. A discussion of the uncertainty 
due to the implementation of the matrix elements with anomalous couplings for the 



multiphoton events is presented in [T0| . 



Table 3: Contributions to the systematic uncertainty on fitted Ak 7 , as explained in the 
text. 



Origin of uncertainty 


Region 1 


Region 2 


Acceptance corrections 


±0.08 


±0.08 


Photon energy calibration ±1% 


±0.10 


±0.20 


Background < 1 event 


+ 0.05 


+ 0.05 


Model uncertainty < ±5% 


±0.10 


±0.15 


Luminosity value ±0.6% 


±0.03 


±0.03 


Total 


±0.20 


±0.30 



The model uncertainty in introducing the anomalous couplings into the simulation 
has been checked. The reliability of the simulation of the Standard Model is 
discussed in [|10|. 



• Another contribution to the uncertainty on the total cross section part of the fit is 
given by the luminosity error. 

Other possible contributions to the systematic error, such as the statistical precision 
on the correction factors for muon rejection and energy deposition in the forward region 
of the detector, are negligible. For A 7 , the basic errors are the same, one region only being 
used for the systematic error calculation. 

7 Conclusions 

The anomalous coupling parameters Ak 7 and A 7 have been measured from single and 
multiphoton events in e + e~ collisions between 161 and 184 GeV. The results from the fit 
to the cross sections and to the energy and angular distributions of the photons are 

Ak 7 = 0.05l^(stat) ± 0.25(syst) 
A 7 = -0.05i{;^(stat)±0.30(syst). 

The corresponding 95% C.L. limits including systematic errors are : 

—2.2 < Ak 7 < 2.3 assuming A 7 = 
—3.1 < A 7 < 3.2 assuming Ak 7 = 

These results are in good agreement with the Standard Model predictions and the 
uncertainty is largely dominated by the limited statistics of the data sample. 
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